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Abstract: We describe a broad class of time-dependent exact wave solutions to 6D 
gauged chiral supergravity with two compact dimensions. These 6D solutions are non- 
trivial warped generalizations of 4D pp-waves and Kundt class solutions and describe 
how a broad class of previously-static compactifications from 6D to 4D (sourced by 
two 3-branes) respond to waves moving along one of the uncompactified directions. 
Because our methods are generally applicable to any higher dimensional supergravity 
they are likely to be of use for finding the supergravity limit of time-dependent solu- 
tions in string theory. The 6D solutions are interesting in their own right, describing 
6D shock waves induced by high energy particles on the branes, and as descriptions of 
the near-brane limit of the transient wavefront arising from a local bubble-nucleation 
event on one of the branes, such as might occur if a tension-changing phase transition 
were to occur. 
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1. Introduction 

Understanding time- dependent dynamics is central to applications of higher-dimensional 
supergravity theories [1] to cosmology and to particle physics. Inasmuch as higher- 
dimensional supergravities provide the low-energy limit of string theories, any under- 
standing of time- dependence in the supergravity limit also provides a guide for the 
thornier issue of understanding these same issues in string theory. For these reasons 
there is considerable interest in finding time-dependent solutions to higher-dimensional 
supergravity [2] (as well as of non-supersymmetric gravity [3], since this can also some- 
times capture similar physics). 

Much of this activity has focussed on 10D and 5D theories (motivated by string ap- 
plications, and Randall-Sundrum [4] constructions), although these can either be more 
difficult to solve or they can have features which are specific to the relative simplicity 
of co-dimension one spaces. Six-dimensional supergravity has more recently emerged 
as being a useful intermediate workshop within which to investigate phenomena which 
can often generalize to still higher-dimensional contexts. Interest in 6D supergravity 
has been further sharpened by the recognition that it can provide insights into the na- 
ture of the cosmological constant problem [5]-[9], by building on the observation that 
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higher- dimensional theories can break the link between the 4D vacuum energy density 
and the curvature of 4D spacetime [10]-[13] (see [14] for a review). There has also 
been considerable recent interest in 6 dimensional models more generally [15]. Includ- 
ing branes is notoriously difficult due the necessity to regularize UV divergences which 
arise [16], however for recent work on understanding this more deeply in the context of 
effective field theory see [17]. 

In this paper we further the program of understanding time-dependent solutions 
to higher-dimensional supergravity in two ways. In §2, we present a method for con- 
structing explicit exact solutions to the supergravity field equations of 6D gauged chiral 
supergravity, to derive a new class of exact solutions to these equations which describe 
a class of gravitational waves passing through spacetimes for which two dimensions 
are compactified in response to the stress-energy of two space-filling 3-branes. In §3 
we discuss the applications of these solutions, including how to construct the shock 
wave metric corresponding to ultra-relativistic particles and BH's moving on one of the 
branes, and also how these solutions provide some insight into the transient part of the 
dynamics describing outgoing waves which would arise shortly after a phase transition 
on one of the source branes. 



2. Wave Solutions to 6D Supergravity 

In this section we describe the technique of finding the wave solutions to the field 
equations of 6D gauged chiral supergravity [18, 19, 1]. 

2.1 Field equations 

The action whose variation gives the field equations of interest is part of the Lagrangian 
density for 6D chiral gauged supergravity, and is given by 1 



2k 2 



I gMN 



Rmn + d M 6 <9 A 



~ e~+ F MN F MN - -L e -*t> G MNP G MNP , (2.1) 

where is the 6D scalar dilaton, G = dB is the field strength for a Kalb-Ramond 
potential and F = dA is the field strength for the gauge potential, A M , whose flux in 
the extra dimensions stabilizes the compactifications. The parameters g and k have 
dimensions of inverse mass and inverse mass-squared, respectively. These expressions 
set some of the bosonic fields of 6D supergravity to zero, as is consistent with the 



1 The curvature conventions used here are those of Weinberg's book [20], and differ from those of 
MTW [21] only by an overall sign in the Ricmann tensor. 
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corresponding field equations. The corresponding field equations are E = E M = E MN = 

E MN = 0, Where: 

E = □ + ^ e 2 ^ G MNP G MNP + ^ e~* E MN E MN -%e* 
6 4 k, 2 

E^ p = D M (e- 2<t >G MNP ^ 

E N = D M (e-* F MAr ) + e" 2 * G MNP F MP (2.2) 

k 2 1 

E MN = R MN + d M cj)d N (j) + — G MPQ G N PQ + ^e^ F MP F N P + - (□ 0) # MJV . 

In what follows we also set G MNP = 0. 

There is an ever-growing literature on the static exact solutions to these equations, 
describing compactifications of 6D down to 4D [18, 5, 6, 22, 23, 24], as well as 4D de 
Sitter solutions [25], time-dependent solutions to the linearized equations [26, 27, 28] 
and exact scaling solutions [29, 30]. 

2.2 Nonlinear pp- waves 

We now consider a remarkable set of exact gravitational- wave solutions to the field 
equations, (2.2). This new class of solutions generalizes the familiar pp-wave config- 
urations to the warped compactifications of the 6D gauged chiral supergravity. More 
precisely they are 6D versions of known 4D solutions belonging to the Kundt class [31], 
for which pp-waves are a special case (see also [32] for other specific examples of these 
types of solutions). There is a large literature on closely related solutions [33]. Such 
pp-waves arise in a wide range of contexts, and have been useful to understand black- 
hole collisions [34], Penrose limits of the AdS-CFT correspondence [35], as toy models 
of cosmological solutions [36], as well as being some of the few known exact solutions 
of the string RG equations valid to all orders in a', i.e. consistent string backgrounds 
[37]. 

For the conventional pp-wave solution the metric takes the form 

ds 2 = 2dudv + H(u,x)du 2 + dx 2 , (2.3) 

where u — x\ + 1 and v — xi — t are light-cone variables, and x collectively denotes the 
remaining spatial directions. The wave profile, H, typically satisfies an equation of the 
form 

V 2 H = J(u), (2.4) 

where J(u) is a source which may contain contributions from w-dependent dilaton gradi- 
ents and form-fields (see [38] for examples of relevance to string theory). Plane waves are 
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a special case of pp- waves for which the wave profile takes the form H = Hij{v)x l x\ rep- 
resenting nonlinear extensions of gravitational perturbations around Minkowksi space- 
time. 

We wish to extend this class of solutions to describe gravitational waves moving 
through the warped compactifications of refs. [5, 6, 22, 23]. To this end we make 
an ansatz motivated by the solutions of ref. [22], but allowing for a more general 
gravitational-wave configuration motivated by the above pp-wave solutions and their 
Kundt class generalizations. For the remainder of this section we adopt ref. [22] 's units, 
for which k 2 — |, and so our ansatz for the metric then is 



where the coordinates (77, 9) parameterize the two internal dimensions, (x±, X2) label the 
directions parallel to the wavefront in the noncompact four dimensions and (u, v) are 
light-cone coordinates along the direction of wave motion. This form is chosen so that 
g vu is the only nonzero component of type g VM . Similarly the only non-zero component 
of type g UM is g uv . This feature considerably simplifies the equations of motion, and 
in particular implies that the inner product, V ■ W, of any two vectors, V M and W N , 
receives no contribution from terms of the form V U W U . Similar properties hold for the 
contractions of tensors of arbitrary rank. 

As it stands this metric is more general than is necessary since coordinate trans- 
formations of the form 



can be used to restrict some of the undetermined functions. In particular we can always 
set N to unity by means of an appropriate rj redefinition, and then set the resulting 
K v to zero by means of a redefinition of v. The net result just gives a redefinition of 
H(r),u,x,9). In what follows we shall set K v to zero, but we do not (yet) fix N, for 
reasons which become clear below. This leads to: 



ds 2 = W 2 {r], u) e ~^' u) 2 du dv + H(rj, u, x, 9) du 2 




(2.5) 



+K ri (ri,u)dr]du + a 2 (r],u) d9 2 + W 8 (ri,u)N 2 (r],u) drj 2 , 





(2.6) 



ds 2 = W 2 {r], u) e -^' u) 2 du dv + H{ V , u, x, 9) du 2 

+W 2 { V , u) e^ u) r e & fo' u > dx\ + e~^ u) dx\ 

+a 2 (r],u)\d9 2 + W 8 (rj,u) N 2 (rj,u) drj 2 . 



(2.7) 
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The fact that H can depend on 9 breaks the axial symmetry corresponding to shifting 
this coordinate, and if a is defined so that 9 goes from to 2rr then it is natural to 
consider periodic functions 



H(r],u,x,9 + 2n) = H(r],u,x,9), (2.8) 

to ensure regularity of the metric at finite t]. We further assume that the dilaton is 
taken to be a function = <f)(i], u) and that the only nonzero component of the gauge 
potential is A e = Ag(rj,u). 

The analysis of the equations of motion is greatly simplified if we make the choice 
of variables [22] 

lna = ^(3X + y + 2Z) 

\nW = ^(y - X) (2.9) 

and so e 2X = aV, e 2y = a 2 W 8 e^ and e~ z = W 2 e*. 

Substituting the ansatz, (2.7), into the field equations, (2.2), we find that the 
nontrivial equations of motion fall into two sets: (i) One set of equations — E, Eg, 
E m , E uv , En, E 22 , and Egg — contain no u derivatives, and so degenerate into ordinary 
differential equations which govern the dependence on 77; (ii) By contrast, the second 
set of equations — E uu and E urj — contain both r] and u derivatives. 

We start with the first set of equations, which involve only derivatives with respect 
to 77. The gauge field equation becomes 

(e-^N-'A'g)' = 0, (Eg) (2.10) 

where ' denotes d/di]. Convenient combinations of the Einstein and dilaton equations 
then give 

(N-'Q'^O {E\-E 2 2 ) 
(N-'Z'Y = (E\ + (E\ + E 2 2 )/2) 
(N- 1 ^)' = (E\ - (E\ + E 2 2 )/2) (2.11) 
N-\N- l X'y + e- 2X K 2 N' 2 (A' e f = (E%) 

N-^N-iyy + %e 2y = (E v - E%/2 - E\ - (E\ + E 2 2 )/2) . 
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Finally the constraint equation becomes 

N- 1 (V 2 - y 2 + Z' 2 + i[ 2 + ^ + « 2 e- 2 *A^ - ^ e 2y iV = 0. (E\ - E M M /2) . 

(2.12) 

This is a constraint equation in the sense that it contains no second-order derivatives 
with respect to 77. In particular, if it is satisfied at any given r\ one can show that the 
other equations imply it must be satisfied for all 77. 

These equations of motion, combined with the constraint, may be derived from the 
following action 



/ 



d?7 



N -l _ y,2 + z ,2 + ^ + ^_ + e ^X K 2 A ^ + ^-e^ N 



. (2.13) 



We now see that our reason for not gauging away the N variable was to use variation 
of N to allow the constraint to be derived from this action. From now on we set N = 1 
in the field equations. 

Of the second set of equations, E ur} is also a constraint equation in the sense that 
it does not contain double rj derivatives. Its explicit form is 

2K 2 e- 2X A' 9 A e , u +2 &i,«+&2,«+2 y'y, u +2 Z'Z, u +^ u +2 y[ u +Z[ u = . (E ur] ) 

(2.14) 

Direct use of the Bianchi identities V M G M N = shows that 

d v E ur) = c AB (r), u)E AB + d(r), u)d u E uv . (2.15) 

Thus provided the other equations are satisfied, and we choose the boundary data so 
that E urj = on a surface i] = r]o, then this equation is satisfied for all 77. 

The final equation (E uu ) is a linear, sourced wave equation for the wave profile H. 
It takes the form 

OH = J(r],u) (E uu ), (2.16) 
where O is proportional to the scalar Laplacian on the spacetime (2.7), 

= d 2 + e~^ +2y+z dl + e ~^ +2y+z d\ + e~ 2X+2y d 2 , (2.17) 

and the source, J, is given by 

J(rj, u) = - e ^ +2y+z [2K 2 e' 2X Al u + 3 £,„ + £, u + 2 X% + 2 y% + Ay, u Z, u 

+3Z% + 2£ 1>u (2y tU + Z tU ) + 2i l)UU + 4y, uu + 2Z tUU j . (2.18) 
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This form is reminiscent of the usual pp-wave equations, except that now the source 
additionally depends on the direction rj. 

The whole system of equations is now straightforward to solve. We first solve 
the ordinary differential equations for (Ag, X, y, Z, £1, £ 2 ), and promote the integration 
constants to arbitrary functions of u. Choosing N — 1, this leads to 

A' e = C^u) e 2X , Z = C 2 (u) V + C 3 (u) 

Si = C A (u)ri + C 5 (u), & = C 6 (u)r ] + C 7 (u), (2.19) 

and 



X" + k 2 C 2 (u) e 2X = 

y + l^ e 2y = o. (2.20) 

K 

The explicit solutions to these equations are 



e y = K^^ S ech(C 9 (u)( V - V2 (u))) . (2.21) 

The restrictions on the integration 'constants', Ci(u), appearing here come from impos- 
ing the constraints — i.e. eqs. (2.12) and (2.14) — on some chosen surface rj = rjo- In 
addition we have the flux quantization constraint which comes from defining the gauge 
field to be smooth on two separate patches and demanding that the associated gauge 
transformation picks up an integer multiple of 2tt phase on integration around 0. In 
practise this is the statement that 

J jF ve d9dr ] = ^, (2.22) 

where q is the smallest unit of fermion charge coupled to A^. This is equivalent to 

Ci(«) e 2X drj = -. (2.23) 



Q 

Substituting in the above solution Eq. (2.21) for X, this relation fixes the function 

Ci(u) to be 

Cl(M) = HiMO. (2 . 24) 

K, z n 

Finally we can compute the source J(i], u) to obtain the wave profile H in the form 

rv rvi 

H(r],u,x 1 ,x 2 ,0) = drji dr) 2 J(t) 2 ,u) + H (r),u,x 1 ,x 2 ,9) , (2.25) 

Jr)o Jr]o 
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where H (r],u,Xi,X2,9) is any solution of the homogeneous equation: OH = 0. Since 
O does not contain u derivatives, in solving this equation we may further allow all the 
integration constants to be arbitrary functions of u. As is the case for normal pp-waves, 
this leaves us free to introduce an arbitrary functional dependence on u into the wave 
profile. 

2.3 Breaking the null symmetry 

A simple interesting extension of the wave solutions can be obtained by allowing for 
the wave profile H to depend on the null coordinate v. Once this is done there is no 
longer a null Killing vector. The ansatz that is consistent with the equations of motion 
in this case is 



ds 2 = W\r],u) e 



2 du dv + H(rj : u, x, 6, v) du* 



+W 2 (r],u)e^ u) 
+a 2 (r],u) 



e 6 (»?,«) d x 2 + e -6(Tj,u) dx 2 
d6 2 + W 8 {r],u)dr] 2 j , 



(2.26) 



with 



H(rj, u, x, 6, v) = if x (77, u, x, 9) + H 2 (r], u)v + H 3 (u)v 2 . (2.27) 

This is identical to the situation for the 4D solutions belonging to the Kundt class 
[31]. From now on we shall focus on the case for which H 3 (u) is independent of u. The 
analysis proceeds as before. In particular, the set of equations involving only derivatives 
with respect to r\ can be encoded in the action 



di] 



N X 



>I2 



y ,2 + z ,2 + ^ + ^ + K 2 e -2X^2 

2 



% 2 



K 



e 2y -2H,e^ +2y+z )N 



(2.28) 



while the E nu constraint becomes 

H' 2 = 2K 2 e- 2X A' e A e , u +2 ^ u+ ^ u+ 2 X' X, u -2 y'y, u +2 Z'Z, u +^ v +2 y[ u +Z[ u , (E +v ) 

(2.29) 

which we regard as an equation to be solved for ff 2 . Substituting this back into the 
E uu equation we find 

OH, = J, (E ++ ) (2.30) 

with the new source given by 

J = -e^ +2y+Z (2n 2 e^Al u + 3 ^ + £, u + 2 + 2 y% + 4 y >u Z, u + 3 Z 2 U 

+2^, u (2y :U + Z, u )+2^ uu + Ay :UU + 2Z, uu ) -d u (t; 1 + 2y + Z)H 2 . (2.31) 

These equations are again straightforward to solve (in principle). 
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2.4 Relation to brane properties 

An important property of the static compactifications is the connection between their 
asymptotic forms near their singularities and the physical properties of the branes 
present at these singularities [23, 39, 25, 29, 43]. These show that in the case where 
gravity is weak the near-brane limit of the extrinsic curvatures are related to the com- 
ponents of the effective 4D brane stress energy. 

Special case: pure tension branes 

In particular, bulk solutions which are sourced by pure-tension branes (in general the 
tension may also be dilaton dependent T = T ((/>)) must have asymptotic limits for 
which the spatial and temporal parts of the metric in the noncompact four dimensions 
are invariant under 4D Lorentz transformations. We next identify the conditions which 
the above solutions must satisfy in order to be sourced in this way by this type of 
Lorentz-invariant brane. 

In this case the symmetry of the boundary conditions requires £i = 0, £ 2 = and 
H' — at the positions of each brane. Inspection of the solutions shows that the first 
two of these conditions require — £' 2 — everywhere throughout the bulk, and so we 
may for convenience everywhere set £i = £2 — 0. On so doing, the remaining equations 
for X, y, Z and Ag are identical to those obtained for the static solutions of ref. [22, 23], 
and we find then the most general solution to have the form: 

£ -(p(r),u) _ yy2 e X 3 (u)i] 



These are precisely the static solutions, but with all integration constants made into 
functions of u, for which the branes are situated at rj = ±00. The choice of 
corresponds to the the fact that we strictly need two patches to define a smooth gauge 
field. As discussed earlier, one function of u (e.g. Q(u)) will be fixed by the flux 
quantization condition which is here the statement that 



W 4 (rj : u) 



( Q{u)X 2 (u) \ cosh[Ai(w)(?7 - 771 (tt))] 
V 4g\i(u) J cosh[X 2 (u)(r] - r] 2 (u))] 



(2.32) 





(2.33) 



and so 



Q(u) 




(2.34) 



n 
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The E m constraint imposes the relation \\{u) = A|(m) + Xf(u). By choosing the 
boundary data for the r\ integration at the surface 77 = ?7 — > — 00 we see that the E vu 
constraint becomes at leading order 

^2^2,^ = A3A3 jM + AiAi )U , (2.35) 

which is just the u derivative of the E m constraint, and is hence automatically satisfied. 
At subleading order we find (Q(u) / Xi(u)) u = which is automatically satisfied if the 
flux quantization condition is satisfied. In other words this constraint does not give us 
any new information on the integration functions. 

The last constraint required for pure-tension branes is the imposition of H' = at 
r] = ±00. In general we see no obstruction to achieving this using the above solution. 
For instance, suppose that H were only to depend on r\ and u. Then we may use the 
homogenous solution to arrange H' — > at one brane. This is consistent with H' — 
at the other brane if the quantity I(u) = 0, where 

/+00 
dr]J(r],u) (2.36) 
-00 

/+00 
d V e 2y+z {2K 2 e- 2X A dtU 2 + 2 X 2 U + 2 y% + 4 y„Z, u + 3Z% + 4 y, uu + 2 Z, au ) . 
-00 

We must first check to see if this integral is convergent. In general convergence of all 
but the first term in brackets is guaranteed by the exponential falloff of e 2y+z at both 
boundaries. On the other hand, we also require A d:U 2 to falls of at the branes as e 2X 
and this is sufficient to guarantee the convergence of the integral as a whole. Given 
that the integral converges, and that the w-dependence of the integration constants 
is arbitrary, it would seem to be possible in general to choose this -u-dependence to 
arrange I(u) = 0, and so H' = 0. 



3. Applications 

The solutions we described in the previous section represent the natural lift of the 4 
dimensional pp-waves and Kundt class solutions up to 6 dimensions, and consequently 
they have many of the same applications as arise for usual pp-waves. 

3.1 Shock waves arising from high energy particles on the brane 

A classic example of a pp-wave is the Aichelberg-Sexl metric [44] which describes the 
metric induced by an ultra-relativistic particle moving through Minkowski spacetime, 

2 

ds 2 = 2dudv - — <f(u) log(|f|)d M 2 + dx 2 . (3.1) 



- 10 - 



It is straightforward to obtain this metric by taking the Schwarzshild metric, which 
describes the gravitational field of a massive particle at rest, and infinitely boosting 
it (i.e. taking the Penrose limit) so that the particle is moving at the speed of light. 
Alternatively we may obtain the metric directly as the pp-wave induced by a stress- 
energy of the form 

T uu = m5(u)5 2 (x), (3.2) 

which is localised on the worldline of a massless particle with energy m. This metric 
is useful in understanding how black holes form in the collision of two high energy 
particles [34]. 

In the present context it is natural to ask, what is the 6D metric describing a 
high energy particle traveling on one of the branes? Not surprisingly the solution is 
given by one of the warped pp-waves we described in the previous section, provided we 
supplement the equations with an additional source term 

T uu = m5(u) 5 2 (x) 5 2 (y-y b ), (3.3) 

where y b denotes the position of one of the branes, and as before m is the energy of the 
ultra-relativistic particle on the brane. The only equation that is modified from before 
is E uu , so that (2.16) is replaced with 

OH = J(r],u) - 2K 2 m5(u) 5 2 (x) 5 2 (y - y b ). (3.4) 

A crucial point is that because this source is traceless (both in a six dimensional = 
and four dimensional sense Tfi — 0), the positions of the branes are not affected by the 
presence of the source. Thus the entire backreaction of the source is encoded in the 
modified wave profile for H. 

Rather than finding the most general solution consistent with this boundary con- 
dition, we shall look for solutions which have an 0(2) symmetry under rotations in the 
x directions, which amounts to choosing £2 = 0, hence Cq(u) = C-j{u) = and are in 
addition 9 independent so that dgH = 0. Both of these are consistent truncations given 
the form of the source. To reiterate the form of the bulk metric is 

ds 2 = W 2 (r], u) e~^ u) [2 du dv + H(rj, u, x, 6) dvfi 

+W 2 (rj, u) e^' u) d 2 x + a 2 (r], u) \d9 2 + W 8 (r], u) drf 



where 



A! e = C 1 (u)e 2X , Z = C 2 (u)r ] + C 3 (u), & = C 4 (u) rj + C 5 (u), (3.5) 
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and 

e x = ^-se C h(C 8 (u)( V - Vl (u))), 

K Ci[U) 

e y = K ^M sech(( C 9 (u)(r 1 -r l2 (u))) ; (3. 6 ) 
^9 

with Ci(u) = 2qC%{u) / K 2 n from flux quantization. We must also check that the 
constraints are satisfied. Equation (2.12) amounts to 

Cl{u) = Cl{u) + Cl{u) + Cl(u), (3.7) 

whereas on substituting this relation into equation (2.14) we get 

'CAu) 



CM 



0. (3.8) 



which is consistent with the flux quantization condition and so does not give us a new 
constraint. The exact nonlinear equation for the wave profile H reduces to 



d 2 -i 1+ 2y + Z^2 



n 



H = J( V ,u) - 2^(# ^-^ , (3.9) 
V/ ' ; V 1 2-nr 2tt ' K ' 



where we write r 2 = x 2 . Because the two branes are located at rj — > ±oo in our 
coordinates, we adopt the artifice of taking the position of the brane to be at i] = 77^ 
with r/b — > ±00, depending on which brane the particle is on, taken at the end of the 
calculation. The expression for the source J(r),u), which is now fully determined, is 
given by (2.18). Let us consider the behavior of the source as rj — > ±00. It is easy to 
see that the leading behaviour is 

J ^ r j2 e (C i {u)+C2(u)) V -2\C 9 (u)\\ v \ ^ 

and so provided |C 9 (u)| > \C±{u) + C 2 (u)\/2, J is bounded at 77 = ±00. Given this we 
can integrate the profile as 

/v rm 
drj! di] 2 J (r) 2 ,u) + H (r],u,x), (3.11) 
-00 J — 00 

where H satisfies 

H = -2^(u)P. W ^r m) . (3-12) 



o 2 v + e -^+ 2 y+^ 2 



2nr 2tt 

The formal solution to this equation can be expressed as the shock wave 

f°° dk 

H (r),r,u) = -2K 2 m8(u) / —kJ (kr)u k (r],u), (3.13) 

Jo 2-k 



- 12 - 



where the modes Uk{f],u) are defined to satisfy 



K 2 C 9 (u) 2 e (C2 {u) ~ Ca ^V+Cs (u)-C 5 (u) 

V cosh 2 (C 9 (u)(r]-r] 2 (u))) 



k 



2vr 



(3.14) 



To see how to proceed let us choose to place the source on the brane at rjb = — oo. We 
then look for a solution to the homogeneous equation which is regular at rj = +00 and 
satisfies the boundary condition 



d Wt 

^Mv,u)\ v=Vb = — . 



(3.15) 



In principal this is straightforward to do numerically. As a specific analytical example 
let us choose C±{u) = 62(14). The general solutions of the homogenous equations can 
be obtained in terms of the variable z = tanh(Cg («)(?? — 772(14))). It is given in terms 
of Legendre polynomials: 



u k (rj,u) = A k P v {z) + B k Q u (z), 



with 



1 

U= 2 



-1 + 



k 2 K 2 C 2 (u) 



,C 3 (u)-C 5 (u) 



9 Z 



(3.16) 



(3.17) 



We can without loss of generality choose Cg(u) > so that the sourced brane is at 
z — — 1. Demanding regularity at z — +1 enforces B k = 0. Note that this is still the 
correct solution even when v is complex, and that P u (z) is still real. The Legendre 
function has the property that 



lim P v (z) 



2 7 + ln(^)+^_ l/ (0)+^ 1+y (0) 



+ 



(3.18) 



with ip the polygamma function. As is usual for codimension two objects, these metrics 
strictly diverge logarithmically at the branes, seen here as the logarithmic divergence 
of the Legendre polynomial at z — — 1. This issue may be understood in terms of the 
regularization and renormalization prescription described in Ref. [17]. For now we shall 
just regularize this by hand. Finally, from the boundary condition we can determined 
A k 

A k = -2C 9 (u)T_ u T l+u -t. (3.19) 

This completes the specification of the solution. 

These metrics could be useful in understanding black hole collisions when the im- 
pact parameter is comparable to the size of the extra dimensions, for which it is neces- 
sary to take properly into account the warping in the extra dimensions. 
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Figure 1: A cartoon of a wave moving into the bulk and along the brane due to a bubble- 
nucleation event on a brane. The equations in the text provide an approximate description 
far from the nucleation site, and near the brane (where the disturbance is approximately a 
plane wave moving parallel to a brane direction). 

3.2 Brane Bubble Walls 

Another physical application of the above solutions is to the problem of understanding 
how 6D compactifications respond to phase transitions on the source branes. It is known 
that such a phase transition can significantly perturb the bulk field configurations, and 
in some circumstances might completely destabilize it. This expectation comes because 
it is known that static solutions can only exist provided that the tensions and other 
bulk couplings of the two source branes are related to one another [5, 40, 23]. More 
recently it has been found that the late-time solutions which arise when the branes are 
not so adjusted are time-dependent scaling solutions [29, 30]. This suggests that if the 
tension on one brane were to be perturbed in an arbitrary way, it could precipitate a 
transition from a static background towards one of the runaway scaling solutions. 

One way to estimate this response is to perform a linearized stability analysis 
about one of the known static compactifications [26, 27, 28]. The response to a change 
in one of the brane tensions may be accounted for in such a calculation by allowing 
the perturbations to become singular at the brane positions in a known way [27]. The 
result of such a calculation is that the perturbations are marginally stable, allowing 
time-dependent deformations to grow along the system's known flat direction. 

This type of linearized calculation is not really satisfactory because it models the 
brane tension as changing simultaneously throughout the brane volume. However 
causality implies that such a transition must occur locally with the nucleation of a 
small region having the new tension at some position on the brane, with the region of 
new tension then expanding to fill space to take more and more of the brane volume to 
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the new vacuum. The resulting geometry is likely to be reasonably complicated, with 
a wave propagating out into the bulk as well as along the brane directions (see Fig. 1). 
We expect the scaling solutions of ref. [29, 30] to be obtained after the wave reflects 
back and forth a sufficient number of times around the internal dimensions. 

The solutions in this paper provide some approximate insight into the nature of 
the waves which would carry the news of the phase transition across the brane and the 
bulk. The insight is only approximate partly because the wave solutions we consider 
are plane waves moving along a direction parallel to the branes, rather than describing 
'spherical' emission into the brane and the bulk from the nucleation event. Our solution 
might nevertheless capture some of the physics of the outgoing wave far from the 
nucleation site, along directions in spacetime near to the brane itself (and so moving 
close to parallel to a brane direction). Since it describes a single wave it also can 
only apply before any waves have had time to be reflected back to the initial brane 
after crossing the extra dimensions. In such a situation we expect the physics of wave 
propagation to be largely insensitive to the specific boundary conditions taken by our 
solutions on the spectator brane. These conditions need not be inconsistent with one 
another, particularly in scenarios for which the extra dimensions are much larger than 
the particle-physics scales which govern the nucleation event itself. 

The change in brane properties with wave passage can then be read off from the in- 
dependence of the integration 'constants' in the pp-wave solutions. For instance, for the 
pure-tension branes the quantities \i(u) appearing in eqs. (2.32) can be related to the 
time-dependent brane stress-energies by formulae such as those given in refs. [39, 29], 
with the result that the 4D tension of the brane located at rj — > ±oo is proportional to 
the quantity 3a± + f3±, where 



4. Conclusions 

Six dimensional supergravity provides a fruitful laboratory for investigating the issues 
which underly higher-dimensional physics in general, and brane approaches to the cos- 
mological constant problem in particular. It does so because 6D is rich enough to 
exhibit many of the properties of still-higher dimensions — like moduli-stabilization 
through fluxes [41], brane back-reaction on internal geometries, chiral fermions and 
Green-Schwarz anomaly cancellation [42] . Yet it is also on the one hand simple enough 
to allow the development of techniques for obtaining physically interesting exact solu- 
tions, but on the other hand not so simple as to be misleading about what happens in 
higher dimensional (in a way which co-dimension one physics sometimes can be). 




and 




(3.20) 
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We have used these properties to explore solution-generation techniques which we 
believe to be applicable to a wide variety of higher-dimensional supergravities. We do 
so by using these techniques to construct a new class of time-dependent exact solutions 
to the field equations of 6D chiral gauged supergravity. These solutions describe the 
physics of nonlinear gravitational waves passing through compactified spacetimes for 
which two dimensions are self-consistently compactified in response to the presence of 
two space-filling branes and a bulk Maxwell flux. 

We believe these methods to merit more detailed exploration. 
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